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Ab stract 

A rapid quasi-three-dimensional analysis has 
been developed for blaoe-to-blade flows in turbo- 
machinery. The analysis solves the unsteady Euler 
or thin-layer Nav ier-Stokes equations in a body- 
fitted coordinate system. It accounts for the 
effects of rotation, radius change, and stream- 
surface thickness. The Baldwin-Lomax eddy- 
viscosity model is used for turbulent flows. The 
equations are solved using a two-stage Runge-Kutta 
scheme made efficient by use of vector uation, a 
variable time-step, and a flux-tased multigrid 
scheme, which are all described. A stability 
analysis is presentee for the two-stage scheme. 
Results tor a flat-plate model problem show the 
applicabil it> of the methoo to axial, radial, and 
rotating geometries. Results for a centrifugal 
impeller anc a rad'al diffuser show that the 
quasi-three-o imen >nal viscous analysis can bt a 
practical design 

Introduction 

Turbomachinery intended to produce large 
amounts of power from a small volume often require 
use of radial-flow or mixeu-flow comoonents, that 
is, components in which the streamwise velocity 
is not strictly axial. Radial-flow turbomachines 
Such as centrifugal impellers, radial diffusers, 
and radi al- inf low turbines have a predominantly 
racial flow direction. Mixed-flow turbomachines 
may oe used wnen restrictions on space prevent a 
completely radial flow. Complicated geometries, 
shock waves, and viscous phenomena make analysis 
of radial- or mixed-flow turbomachines more dif- 
ficult than analysis of strictly axial-flow 
machines . 

Analysis of axial-flow turbomachinery blade 
rows is usually simplified by modelling a blade 
section as a flat cascade. The governing equa- 
tions for a flat cascade are the same two- 
dimensional flow equations that are solved for 
isolated airfoils, so cascade analyses often jra* 
heavily on numerical techniques developed for 
two-dimensional external flows. Examples of flat 
cascade analyse include panel methods, poten- 
tial methods , 1 Euler solutions, anc Navier- 
Stoxes solutions .5,6 

Flows in radial- or mixed-flow turbomachines 
are inherently three-dimensional , requiring speci- 
fication of the axial, radial, ana tangential 
velocity components to fully specify the flow. A 
simplification that allows these machines to be 
analyzed in two di.nensions mas proposed by 
C.H. Mu in 1952. In Wu's model the flow is 
assumed to follow an axisymmetric stream surface 
(Wu's M S2 surface," Fig. 1). The radius and 
thickness of the stream surface are assumed to be 
know as functions of the streamwise distance. 

These quantities are usually obtained from an 
ax isymmetr ic through-flow or "meridional* 


analysis,® sometimes coupled with a boundary- 
layer analysis 9 on the hub and shroud. 

The equations governing the flow along the 
stream surface combine the axial- and radial- 
velocity components into one streamwise component, 
and are thus two-dimensional. The solution can be 
resolved into three velocity components since the 
shape of the surface is known. Specifying the 
stream-surface thickness aUows variable Place 
heights and end wall displacement thicknesses to 
bemooclled. This is similar to specifying area 
change in the one-dimensional nozzle equations. 

Since the effects of raoius change and stream- 
surface thickness are modelled in mis analysis, 
it is termed "quesi-tnree-oimensional . " Examples 
of quasi -three-aimens ipoal turbomachinery analyses 
include panel methods, stream-f unct ion 
methccs,* 1 potential methods,** ana Euler methods. ^ 

In the present work, the Euler and Navier- 
Stokes code developed for flat cascades in Ref. 6 
has been extended to a quasi-tnree-oimensional 
analysis. It is thought that this is the first 
Nav ier-5tokes analysis to include the effects of 
rotation, radius change, and stream-surf ace thick- 
ness. The explicit MacCormack algorithm^ used 
in 6 has been replaced with an explicit two-stage 
Runge-Kutta f inite-d inference algorithm based on 
the work of Jamosc^ ** Efficiency is achieved 
by three means; vec*orizetict , use of a variable 
time-step, ^nd by use of a multigrid scheme devel- 
oped bv Ni*- and modified by Johnson and 

Chima .*7-19 

Governin g Equations 

The axi s/rmetric (m,e ) cooroinate system 
uS'd for the quasi-three-oimensional analysis is 
shown in Fig. ], Here the m-coordinate is 
aef in *0 by 


dm * qz ♦ or v 1 

and the o-coordinate is defined by: 

• ■ - -t (2) 

where a' is fixed in space and e rotates 
with the blaoe row wit'i angular velocity r*. The 
radius r and the stream surface thickness h 
are taken to be known functions of m. In this 
coordinate system the dimensionless Navier-Stokes 
equations may be written in the following nearly- 
conservative form: 

» t Q ♦ » m (F - «e' l R) ♦ > e (G - Re' l S) -PC (3) 
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*6 * v e “ ^ ** the relative tangential velocity 

component, 

e - p[CyT ♦ l/2(v£ + vf)] is the total energy 
per unit volume, 

p « U - 1 )[e - 1 / 2 ( vf* ♦ vf )] is the pressure, 

* (ev^ ♦ p - o^Jrm/r ♦ (p _ a 33 }hm/h 


r /r ■ - — 
m r dm 


n /h . i 
m n dr 


The viscous terms in the energy equation are: 
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a = Y>P/c is the sonic velocity, and the nor- 
malised thermal conductivity k equals one. 

Tne shear stress terms are: 
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The equations are nonaimensional i 2 ed by arbi- 
trary reference quantities (here the inlet total 
density and critical sonic velocity define the 
reference state), ana the Reynolds number Re and 
the Pranatl number Pr must be specified in terms 
of tnat state. These equations assume that the 
specific neats C p and C v ana the Pranctl number 
are constant, that Stokes' hypothesis x . -2/3 u 
is valid, ana that the effective viscosity may be 
written 


u lamirar * ^turbulent 


Equations i3) to (6) are transformed from the 
(m,a) coordinate system to a general bocy-fitted 
U,n) coordinate system using standard methods. 

Tne thin-layer assumption is then used to elimi- 
nate viscous aerivitives in the streamwise (c ) 
direction, thereby reducing computational overhead 
while retaining the capability cf computing sepa- 
rated flows. The resulting equations are similar 
to those developed by Katsanis.™ 
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Computational Grid 


K 2 " ( BV * 2 * p ' Re ' l °22) r n/ r * <P - Re ' 1 °33 ) V h 

s « - Tt - k n*»r (% ♦ % )v 2 * V? * V3 

In Eq. (8) the overbar: denote a rescaling 
of the metric terms: 


** “ *» /r; % ■ V r; J_1 • rhJ * 1 (*) 

where J is the Jacobian of the transformation 

J - *■% - v. ■ r ^> r> - i^T (10) 

and the metric quantities are determined from the 
grid-point coordinates using central clfferences 
and: 



The relative contravarlant velocity components 
* and U along the c and n grid lines are 
given by: 
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The shear stress terms are found from Eq. (6) 
by replacing a m with n m i n and l/r* 6 with 


The quasi-three-dimenslonal equations 
(Eqs. (7) to (12) ) are similar to the two- 
dimensional equations solved in Ref. 6 except for 
the source term Kj* the radius appearing in the 
•-momentum equation, the rescaled metrics 
Uq. (9)), and the relative velocity component 
W| appearing in the contravarlant velocities 
(tq. (12)). Equations (7) to (12) reduce to the 
two-dimensional equations for constant r and h, 
and zero rotation. Note that Eqs. (7) to (12) are 
independent of the magnitude of the stream surface 
thickness h so that any function h(m) > 0 may 
be used. The equations do depend on the magnitude 
of the radius r because of the 1/r terms scal- 
ing c e and n e in Eq. (9). 


Body-fitted grids for this work were 
generated using the GRAPE code (GRids about Air- 
f Jils using Poisson 1 * Equation) developed by 
Sorenson/**^ Briefly, the code allows arbi- 
trary specification of Inner and outer boundary 
points, then generates interior points as a solu- 
tion to a Poisson equation. Forcing terms In the 
Poisson equation are chosen such that desired grid 
spacing ana Intersection angles may be maintained 
*t the inner and outer boundaries. 

hew Inner and outer boundary subroutines were 
written for turbomachine geometries. The new 
inner bounoary has constant spacing around blade 
leading and trailing edges, larger'constant spac- 
ing over blade surfaces, and exponential stretch- 
ing connecting the regions. The new outer 
boundary Is composed of the mean-camber line 
between the blades, a quadratic extension upstream, 
and a linear extension downstream. The C-snaoed 
grids are periodic over the pitch of the blade. 
Grids are generated in a Cartesian (m,7s) coor- 
dinate system, where 7 is some mean *ad1us. 

The local radius and stream surface t «iickness are 
supplied to the Navier-Stokes :^ue at a later time 
as tabulated functions of a, then spline-fit on 
to the grid. The * . r m /r and h^/h In 
Eqs. (7) and (8) are calculated using central 
differences and are stored. 

Initial Conditions 


Since a centrifugal compressor can produce 
pressure ratios of 5:1 or greater it Is not gen- 
erally possible to start a quasi-three-dimenslonal 
calculation with constant Initial conditions. 
Instead an analytic solution of the one-dimensional 
flow equations with area change is used. Turoo- 
machlnery blades are usually designed based on 
desired leading- and trai l ing-edge velocity tri- 
angles, so the relative flow velocity w and 
angle a at the leading edge and the relative 
flow angle at the trailing edge are used as Input 
for the Initial conditions. 

Using ()' to denote absolute total condi- 
tions and ()" to denote relative total condi- 
tions, the continuity ano energy equations are 

m ■ prhaekl cos a * constant (15) 

I « CpT' - n,v ft ■ C p T w - 1/2 r 2 fl 2 « constant 

(16) 


For turbulent flows the two-layer eddy- 
viscosity model developed by Baldwin and Lomax^l 
is used. In the (m,e ) coordinate system the wall 
shear t* and vorticity » required by the model 
are given by 

’ °12 - ‘ " (V. + 1/r Vm * Vm'') w < 13) 
“ * 7 (v. - 1/r Vm + v .V r ) no 


where ae >s the blade spacing and I is the 
rothalpy. 

Outside of the Made row the angular momentum 
equation gives: 

rv e . constant ( 17 ) 

Equations (16) and (17) can be used with the isen- 
tropic relations 1.0 show that total conditions are 
corstant outside of the blade row. Evaluating 
Eq. (15) at some point () 2 and using £q. (17) 
to eliminate a gives: 
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Now the isentropic relations and the defini- 
tion of t 2 m can be used to eliminate p 2 giving: 


The axi symmetric m-momentum equation may be 
written as* 
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Substituting Eq. (19) into Eq. (18) gives: 
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is the upstream-running Riemann invariant. For 
Steady axisymmetric flow Eq. (22) may be written 




rjr ♦ av h /h 
m m m 


Equation (23) is backward-differenced and 
solved for R* at the inlet. Then v m is found 
from R and the specified whirl rv e and total 
temperature T‘ using: 9 


(v-l)R~ * V( T *l)(4C n T- - 2v 2 ) - 2(v-l)(R-) 2 


t>2 r 2^2 

Equation (20) is solved for w ? at each 
grid point upstream of the blade usfng Newton 

usfna Vo ^ qu4ntit1es then found 

tnp (1 ^* tbe know '' tot4) conditions, and 

tne lsentropic relations. 

F „ “iV 1 "® th ® bl4de ro " rv 9 * constant, so 
EQ- (12) is replaced with an assumption that the 

HZ* 1 a S 9 l e /* 2 , V4ries linearly through the blade 
row. A derivation similar to that above nives- 
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which is solved at each grid point within the 

blade row. Once the flew conditions are known at 

rtn C n*» ai ,n9 e ^ 9a ’ Eq - (20) C4n be used for the 
downstream region. 

Bo undary Conditions 

At the inlet, total pressure, total tempera- 

I* 9 are s P ecif 1ed - For suosonic 
inflow the governing equations have one negative 
Eigenvalue so that one variable at the inlet must 

actprTfr 6d a ? f art ° f the solut,on - Here a chal 
acteristic relation is used to extrapolate the 

upstream-runninq Riemann invariant to the inlet. 


Density and energy are founo using isentropic 
relat ons. 

t I'f ex j, t the st4tic Pressure is specified 
and the other flow quantities are found using 
first-order extrapolation. Zeroth-order extrapo- 
' 4b, ° n ** not usually sufficient in cylindrical 
coordinate systems because the radius and thus 
the velocity may change between grid points. For 
the same reason the inlet and exit boundaries 
cannot be placed arbitrarily far from the blades. 
Too small a radius can cause t'.e flow to be super- 
sonic and to large a radius ca 1 caus* the velocity 
to approach zero at the bounder fes. y 

Blade su*r'ace pressures are found from the 
normal momentum equation: 

U -"m* VeV + ^ n m + 9 e'* n b 


+ •■•‘V. - Vm> V r 


where W 4 . 0 on the surface for viscous flows. 

Periodic boundaries are solved like interior 
pu 1 n l ) « 

Fine-Grid Algorithm 

On the fine grid an explicit two-stage 

[he’!n;? 1 !f e, ; enCe R V2 q ?- Kutt4 algorithm based on 
the work of Jameson 15 is used. It is qlven 
below as applied to Eq. (7). 
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Two-stage Runge-Kutta algorithm 

*(1 ) *n if0 n 
q * q a^AtR 


where 


1.2 


S l2 > .3" - 


(26) 


* - Be ' 1 
Arti f i c i a 1 dissipation 

q n+1 . $ (2) * 0 ‘ 2) 

Collect residuals 

*n + l *n+l *n 


(27) 


aS" 1 1 ft 

Aq = q - q 


(28) 


The two-stage scheme given by Eq. (26) has a 
Courant number limit of one. It is used in pref- 
erence to a hiqher-order scheme with a higher 
Courant number limit because the multigrid scheme 
used here also has a Courant number limit of one. 


Four minor advantages of the Runge-Kutta 
scheme over the MacCormack scheme are noted: 

1. A steady Runge-Kutta solution is independ- 
ent ot the time step while a steady MacCormack 
solution is not. This is not true when the arti- 
ficial dissipation is added in a fractional step 
as in Eq. (27), 

2. The central ly-di f f erenced Runge-Kutta 
scheme produces perfectly synmetric solutions for 
symmetric problems while the one-sided MacCormack 
scheme does not. 


3. for a two-dimensional problem with 
centrally-differenced metrics, the Runge-Kutta 
scheme exactly conserves a free stream while the 
MacCormack scheme does not. Although the Runqe- 
Kutta scheme is not fully conservative for the 
quasi-three-dimcnsional problem because of the 
source term, it has been found to possess better 
conservation properties in general. 

4. The Runge-Kutta scheme is slightly easier 
to program than the MacCormack scheme. 

Differences in convergence rates between the two 
schemes are negligible for Courant numbers near 
one. 


Artificial Dissipation 


Dissipative terms consisting of second and 
fourth differences are added to prevent odd-even 
point decoupling and to How shock capturing, 
the dissipative terms are similar to those used 
by Jameson^ and others. A one-dimensional 
version (c-direction) is given below. In two 
dimensions the dissipation is applied as a 
sequence of one-dimensional operators. 


where 


C 



-2 


v 2 - 0(1) 

U4 - 0(1/16) 

The terms in the coefficient C balance 
similar terms in Eq. (26). In smooth regions of 
the flow the dissipative terms are of third order 
and tnus do not detract from the formal second- 
order accuracy of the fine-grid scheme. In 
regions of the flow where the second difference 
of the pressure is large, the second-difference 
dissipation becomes locally of first order. Note 
that in other work including Ref. 15 the term 
iPrri is commonly divided by an average pressure. 
This is not done here because pressures through a 
centrifugal compressor can increase by factors on 
the order of five, which would decrease the 
dissipation correspondingly. 


jtabi 1 i ty Analysis 

A stability analysis of the fine-grid algo- 
rithm is performed in two parts. The first part 
examines the mooel problem considered by Jameson^ 
and is used here to choose the parameters in the 
two-stage scheme. The model problem is the one- 
dimensional convection equation with third-order 
artificial dissipation: 

<U + <lx * uA* 3f \ xxx « 0 (30) 


Applying the two stage scheme (Eqs. (26) and (27)) 
gives: 

q (1) ■ q" - »! 7 vq j+1 - q..j) n 

q<2) -«•" "7<Vl - q M )(1) 

q n + 1 - q (2) - (V2- 4q j-l +6 V 4q j + l*V2 ) (31) 

where a = At/Ax is the Courant number. 


If we consider a Fourier component of the 
solution 

^ = vV px . vV* (32) 

wnere V is the amplitude, p is the wave number, 
K « pax, and i * y^T, then the amplification 
factor is given by: 

, 

9 « ~r ' ll - 4,1(1 - cos tn 
V 

(33) 

2 2 

x (1 - Ai sin i - a ^ a sin t, ) 

= damping factor x characteristic 
polynomial for undamped scheme. 
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The damping coefficient » is chosen such that 
g - 0 at t w «, giving u - a). 

An undamped n-stage scheme can only be 
staDle to a • n - 1, so x Is taken to be one 
while is chosen. Figure 2 shows a plot of 

Eq. (33) for several values of with x • 1 
and u ■ 1/(16 a). It can be seen that the two- 
stage scheme is stable for 0.5 < *\ £ 1.7, with 
«i « 1.2 giving the best damping over the range 
or frequencies. For * 1.2 and » ■ 1/(16 a) 
it can be shown that the two-stage scheme is 
stable for Courant numbers a & 1.1. In general 
the two-stage scheme is first-o~der accurate in 
time. It is second-order accurate in time only 
if a} * 0.5 . 


the coarse grids. Thus the multigrid scheme used 
here is based solely on the Euler equation!. It 
is entirely independent of the viscous terms, the 
turb'i’ence model, and the artificial dissipation 
us<*d on the fine grid. 

One-step lax-Wendroff sc' nes including N1*s 
scheme require temporal derivatives of the flux 
vectors. These terms are computed as the Jacobian 
matrix of the flux vector times the temporal dif- 
ference of the solution vector. Johnson replaced 
these lengthy computations with a direct temporal 
difference of the flux vector using the old and 
new solutions on the fine grid. 18 This "flux 
based" scheme is considerably simpler than Mi's 
original scheme. 


The second part of the stability analysis 
consider the linearized Euler subset of the 
governing equations (Eqs. (7) and (8)). A 
Von Neumann analysis shows the stability limit on 
the lime step to be: 


The flux-based multigrld scheme is derived 
by expanding the fine-grid change aq n , 

(Eq. (28)), :n a Taylor series. 

aq 0 * 2 » &q n41 ♦ at(aq 1+ ^) t ♦ 0(at 2 ) 


at < 


* + + a V U * * + X ^ 


vjJr 


2a V u * + n " ) " * u > * ^ 

which is implemented as: 


-i 


Tne Euler equations are used to reolace the third 
term. 


(4q" +1 ) t - [at(K - F t - G n )] t (35) 

Interchanging the space and time derivatives and 
using bacxwaro differencing in ti me gives: 


1 v J 
m 


at < CFL x I i vj dm ♦ iw^io* ♦ a 


where 


VI 


♦ —~== 

2 #y<lm 2 ♦ 


dm a i e i ♦ i n • 
s m m 


(34) 


<» Z • 'V + ' V 

A spatial ly-variable time-step At, j is 
used in the two-stage scheme (Eqs. (26) ind (27)) 
to accelerate convergence to steady state. 
Equation (34) is used to choose the time step at 
each grid point such that the Courant number is 
constant, typically CFL - 0.95. Time steps are 
calculated based on the initial conditions. They 
are stored and not updated During the 
calculations. 


Multigrid Algorithm 

The multigrid algorithm originated by Ni 18 
and modified by Johnson and Chima 17 " 19 is used to 
accelerate convergence of the fine-grid algorithm. 
Ni's scheme is basically a one-step Lax-Wendroff 
schema applied on a coarse grid. N1 used it to 
accelerate his own fine-grid Euler scheme. 8 
Johnson adapted N1*s method to other f Joe-grid 
schemes including MacCormack * s scheme. 17 He 
also used it for viscous flows by demonstrating 
that dissipative terms need not be Included on 


<**"* l >t - *t ( K t - r u - t u ) 

. (K ntl - K n MF n+1 - F") r (G " +1 - G") n 
and finally 

aq n+Z - 43" n ♦ at[a* - UF) { - u£) n ] n+1 (36) 

where 

4 K n+1 - K(q' 1+1 ) - K ( q n ) , etc. 

Equation (36) is implemented on a coarse grid 
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spacing 

i 

, and 

ian. 

and time step 

‘ 4t i. 
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* • * 
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4t i.j 
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(37 


Fine-grid changes 

aq n+1 are 

restricted to 


the coarse orid by Injection. The multigrid algo- 
rithm (Eq. 137 ) ) computes coarse-grid corrections 
a!} n * 2 based on changes in the fine-grid solution. 
Thus the coarse-grid corrections must vanish if 
the fine-grid solution converges, thereby retain- 
ing fine-grid accuracy. 
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The coarse-grid corrections are prolonged 
back to the fine grid using bilinear interpolation 
and the fine-grid solution Is updated. The 
process may then be repeated on a coarser grid. 

Vectorlzation 

The explicit Runge-Kutta and multigrid algo- 
rithms used here have been highly vectorized for 
the Cray I-S at NASA Lewis Research Center. 

Indeed the Runge-Kutta computations were clocked 
at about 40 million floating point operations/sec 
(40 mflops) for an Euler solution on a 113 by 25 
grid. The efficiency of the multigrld computa- 
tions decreases as the grid gets coarser and the 
vectors get shorter. 

Tne code was redimensioned for each grid size 
run and required 260 K words of memory for the 
largest grid (161 by 33). The quasi-three- 
aimensional code requires about 20 percent m or* 
storage than the two-dimensional code. 6 

Results 

Results are presented for the following 
problems: a cascade of thin flat plates with 

round leading edges, a centrifugal Impeller, and 
a vaned-radial diffuser. Both Euler and havier- 
Stokes results are presented. 

To aid in developing the quasi-three- 
dimensional code and to illustrate the capabili- 
ties of the analysis, a model problem was devel- 
oped representing a cascade of thin flat plates 
with round leading edges. The plate has unit 
chord, four percent thickness, and a pitch of 0.7. 
Figure 3 shows the computational grid around the 
plate. The inviscid flow grid had 113 by 25 
points with 5* spacing around the leading edge 
circle and a normal spacing of i.8xl0 . The 
viscous flow grid had 113 by 33 points with a 
normal spacing of 5xl(T 5 . The grid shown In 
Fig. 3, and most of the subsequent grid and con- 
tour plots in this paper, are drawn in a t ans- 
formed plane in which the abscissa is / dm /r 
and the ordinate is a. This transformation 
preserves angles and is discussed in Ref. 10. 

Grids are generated in an (m,re ) coordinate 
system where r is some mean radius, and are 
independent of the local radius r and stream- 
surface thickness h. Values of r and h are 
supplied to the quasi-three-dimensional code later 
and can be varied to simulate different geometries. 

A flat cascade was simulated by setting 
r « constant and h * constant. Figure 4 shows 
static pressure contours for an Euler solution 
with an inlet Mach number of 0.33. (Titles on 
the figures give nominal value r . of Mach number, 
Reynolds number, and flow angle based on the 
leading-edge velocity triangles g*ven as input. 
Actual inlet conditions may change as the solution 
develops.) The contours show the stagnation 
region and a mild acceleration due to blockage. 
Identical results were obtained using the two- 
dimensional code. 6 The solution was run 5000 
cycles with two multigrid levels and took 165 sec 
on the Cray. The residuals were reduced eight 
decades so the solution could have been stopped at 
about 2000 cycles with a three-decade reduction 
in the residuals. 


A radial duct with constant crcss-sectional 
area was simulated by setting r * m and 
r x h • constant. The solution is identical to 
the flat cascade results shown in Fig. 4. 

A radial diffuser with 36 blades was simu- 
lated by setting r - m and h * constant. 

Figure 5 shows static pressure contours for an 
Euler solution with an Inlet Mach n*nnber of 0.38 
and an exit Mach number of 0.27. The plot is 
shown In polar coordinates. The contours show the 
diffusion of the flow due to the area change. The 
solution was run 2500 cycles with two multigrid 
levels and took 92 sec on the Cray. The residuals 
were reduced 6 decades so the solution could have 
been stopped at about 1250 cycles with a three- 
decade reduction in the residuals. 

A centrifugal impeller with 36 blades was 
simulated by setting r * m, h * constant, and a 
rotation rate of 1000 raa/sec (9550 rpm). 

Figure 6 shows Mach number contours for a tur- 
oulent solution with an inlet Mach number of 0.5 
and an exit Mach number of 0.41. The impeller 
produces a total pressure ratio of 1.05. The 
plot is shown in polar coordinates with rotation 
upwards. The inlet whirl is zero but because of 
the rotation the in et relative flow angle is 
10.6 . This incidence desymmetr izes the flow ano 
produces e pressure loading on the blade. At the 
exit this loading conflicts with the specified 
exit pressure and Is responsible for the kinks in 
the contours there. 

Figure 7 shows convergence histories for the 
previous example with various levels of multigrid. 
Convergence is taken to be a three decade drop In 
the maximum residual to 2x10 . For this example 

two grids are 1.64 times faster than the fine grid 
and three grids are 3.12 times faster than the 
fine grid. CPU tlmss are Included on the figure. 

The remaining results are for a nominally 6:1 
total-pressure ratio centrifugal Impeller and a 
matching vanea-radial diffuser. These components 
were designed by the Air Force for use in an 
auxilliary power unit and were modified and tested 
at NASA Lewis Research Center. The compressor 
has a diameter of 16.1 cm. It has 19 blades and 
was designed for a total pressure ratio of 5.9:1 
with a mass flow of 1.033 kg/sec at 68 384 rpm. 

The diffuser has 27 blades and a diameter of 
25.1 cm. Further details concerning tne compon- 
ents can De be found in Ref. 25. 

Figure 8 shows the computational grid used for 
the impeller. The grid has 161 by 33 points (161 
by 17 shown) with a normal spacing of 3xlQ~* cm 
for inviscid flows and 5xlCr 6 cm for viscous 
flows. The leading edge is round with a radius of 
0.045 cm and the grid points are 7* apart. 

Figure 9 shows normalized radius (RMSP * r), 
stream-surface thickness (8ESP - h), and their 
product (R*B * rh) versus fraction of impeller 
chord. These values were taken from a meridional 
analysis. 8 

Euler and Navier-Stokes solutions were each 
run 2000 cycles with two multigrid levels, reduc- 
ing the maximum residual three decades. The Euler 
solution took 117 sec ana the Navier-Stokes 
solution took 185 sec on the Cray. 
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Surface static pressure distributions for the 
Impeller are compared In Fig, 10* Three solutions 
with Identical mass flows are shown: an Euler 

solution (dashed), a Navier-Stokes solution 
(solid), and a panel solution (circles, 10 ). The 
Euler and Navier-Stokes solutions were run with 
various exit pressures until the deslrtd mass flow 
was obtained. The shapes of the three pressure 
distributions are similar but the panel solution 
has higher pressure levels since It is loss-free. 
The Euler solution has strong shock losses which 
lower the pressure levels. The Navier-Stokes 
solution has weaker shocks due to viscous smooth- 
ing of the leading edge, but blockage eff' :ts 
decrease the pressure levels overall. 

Figure 11 shows relative Mach number contours 
for the Euler solution. The dashed line Is the 
sonic line. At the inlet the relative Mach number 
is 0.91. At the leading edge the flow has 10 # of 
incidence. This produces a large supersonic 
bubble with a peak Mach number of 1.92 on the 
suction (upper) surface. The bubble terminates 
with a normal shock that is smeared due to grid 
shearing in this region. There is also a tiny 
supersonic bubble on the pressure (lower) surface 
which is not visible at this scale. 

Some interesting flow phenomena evident in 
Fig. 11 can be explained by the concept of a 
relative eddy. The flow through an impeller Is 
predominant ly inviscid and tends to remain irrota- 
tional. The blade row in Fig. 11 is rotating 
downward and so adds clockwise vorticlty to the 
flow. To remain irrotational the flow develops a 
counterclockwise circulation within the passage. 
Thus the flow can be modelled as a superposition 
of a through-flow component and a component rotat- 
ing counter to the blade row called a relat/e 
eddy. The effect of the relative eddy is to 
accelerate the flow on the suction surface and 
decelerate the flow on the pressure surface as 
can oe seen in Fig. 11. 

A more dramatic effect of the relative eday 
is to sweep the unconfined flow beyona the trail- 
ing eoga up and away from the blade in a spiral, 
with a slip line leaving the trailing edge. As a 
common example of slip from rotating machinery, 
consider the flow of wattr from a lawn sprinkler 
as viewed rotating with the sprinkler. The /low 
spirals up and opposite to the rotation of che 
sprinkler. It is emphasised that this is strictly 
an inviscid phenomena. 

Figure 12 shows relative Mach number contours 
for a Navier-Stokes solution for the impeller. 

Here viscous effects reduce the peak suction- 
surface Mach number to 1.72 and the peak pressure- 
surface Mach number to Just over 1.0. Both 
surfaces have small leading-edge separation bub- 
bles that are barely resolved on this grid. The 
pressure surface quickly develops a thick boundary 
layer and the suction-surface boundary layer 
thickens after the shock. This blockage causes 
the viscous pressure levels to be lower than the 
inviscid levels in Fig. 10. The rotation ener- 
gizes and diminishes the boundary layers on the 
radial portions of the blade. 

Here again the relative eody sweeps the flow 
off the trailing edge in a spiral. Since the flow 
does not follow the grid lines the thin-layer 
assumption may be Invalid and diffusion across tfe 


wake may not be properly accounted for. This is 
one shortcoming of the present analysis. However, 
since the Lraillng-edge slip is an inviscid 
phenor ina it Is felt that the character of the 
solution Is correct. 

Figure 13 shows the computational grid for 
the radial diffuser vane. The grid has 145 by 33 
points (145 by 17 shown) with a n:-ma1 spacing of 
2.5x10“* cm for Inviscid flows and 6x10“^ cm for 
viscous flows. The round leading edge has a 
radius of 0.025 cm and the grid points are 7* 
apart. At the trailing edge the actual vane is 
cut off at constant radius and the duct turns 
axially. For this analysis the trailing edge was 
sharpened and the duct was extended radially. 

Figure 14 show:* the normal izea radius ar.d 
stream-surf ace thickness versus fraction of vane 
chord. These values were taken from a combined 
meridional /boundary- layer analysis. 8 * 9 Although 
the vanes have constant height, boundary layer 
blockage decreases the flow area by nearly 
50 percent. 

An Euler solution was run 4000 cycles with 
two multigrid levels, taking 214 sec on the Cray. 

A Navier-Stokes solution was run 2000 cycles with 
three multigrid levels, taking 179 sec. In each 
case the maximum residual was reduced three 
decades. 

Surface static pressure distributions for the 
diffuser vane are compared in Fig. 15. Again, 
panel, Euler, and Navier-Stokes solutions are com- 
pared. The shapes of the pressure distributions 
are similar but the panel solution has higher 
pressure levels since it is loss-free. The Euler 
and Navier-Stokes solutions each have small super- 
sonic bubbles terminated by normal shocks near the 
leading edge, and the losses lower the pressure 
levels. Additional blockage effects cause the 
viscous pressure levels to be even lower than the 
inviscid levels. 

Figure 16 shows Mach number contours for the 
Euler solution. The inlet Mach number is exactly 
1.0 but the radial component is only 0.29 so the 
flow is subsonic in character. Because of the 
increasing radius the total Mach number drops to 
about 0.83 near the leading edge. The flow has a 
slight positive incidence at the leading edge and 
the flow accelerates to a peak Mach number of about 
1.5 on the suction surfore. The tiny supersonic 
bubble (visible only as a black dot at this scale) 
terminates with a shock at about 0.26 percent 
chord. The pressure surface develops an even 
smaller supersonic bubble. The diffusion through 
tne passage is evident in Fig. 16. 

Figure 17 shows Mach number contours fot the 
Navier-Stokes solution. Here viscous effects at 
the leading edge entire^ suppress the formation 
of supersonic bubbles. The thick boundary layers 
that develop in the diffuser are obvious, but the 
flew remains attached on both surfaces. 

Con ducing Remarks 

A quasi-three-oimensional Euler ana Navier- 
Stokes anal vis technique has been developed for 
D lade-to-b lade flows in turbomachinery. The 
analysis solves tne nhin-layer Navier-Stokes 
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equations written ir. general coordinates for an 
axi symmetric stream surface, and accounts for the 
effects of blade-row rotation, radius change, and 
stream- Surface thickness. It is believed that 
this is the first Navier-Stokes analysis to 
include these effects. 

The solution technique Is a two-stage Runge- 
Kutta scheme based cn the work of Jameson, Effi- 
ciency is achieved through use of vectorlzation, 
a spatial ly-variable time-step, and a multlgrid 
scheme based on Johnson's revisions of Ni's 
scheme. Thj multigrid scheme typically reduces 
the CPU time required by Che fine grid scheme 
alone by a factor of about three, for both 
inviscid and viscous mows. 

Results for a model problem show the analysis 
to be viable for a variety of axial, radial, and 
rotating geometries. Results for a centrifugal 
impeller and a radial diffuser vane show that the 
analysis can predict a number of phenomena that 
are not accounted for in previous analyses. These 
phenomena include: leading-edge stagnation 

points, leading-edge separation, supersonic 
regions and shocks, blade-surface boundary layer 
growth, and trai 1 i ng-edge slip lines. 


7. Wu, C.-H,, "A General Theory of Three- 
Dimensional Flow in Subsonic and Supersonic 
Turbonachines of Axial-, Radial-, and Mixed- 
Flow Types," NACA TN-2604, 1952. 

8. Katsanis , T., and McNally, W.D., "Revised 
FORTRAN Program for Calculating Velocities 

* nd Streamlines on the Hub-Shroud Midchannel 
Stream Surface of an Axial-, Radial-, or 
Mixed-Flow Turbomachine or Annular Duct, I - 
User's Manual," NASA TN-D-8430, 1977. 

9. McNally, W.D*, "FOR T RAN Program for Calculat 
ing Compressible Laminar and TurbuTent Bound- 
ary Layers in Arbitrary Pressure Gradients," 
NASA TN-D- ; C61, 1970. 

10. McFarland, E.R., "A Rapid Blade-to-Blade 
Solution for Use in Turbomachinery Design, M 
Journal of Engineering for Gas Turbines and 
^owrr, ToTT 1.06. No. j. Ap rTljgTT 

pp. T76-382. 

11. Katsanis, T. , "FORTRAN Program Ter Calculat- 
ing Transonic Velocities on a l lade-to-B lade 
Stream Surface of a Turbomachi ne," NASA 
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It is thought that the ability to predict 
these phenomena rapidly for general geometries 
could make the quasi-three-dimensional analysis a 
useful tool for turhomachinery design. Further- 
more the quasi-three-dimensional analysis can 
provide insight into both physical and numerical 
problems that can be expected with fully three- 
aimensional problems in the future. 
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Figure 1. - Ouasi-three-dlmer sional stream surf,*ce and coordinate 
system for a centrifugal compressor. 



Figure 2. - Ampliiication factor for two-stage Runge-Kutta scheme. 
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Figure 3. - Computational grid (or flat plate model problem. 
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Figure 4. - Static pressure contours for Invisckl flat plate cascade model 
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Figure 5. - Static pressure contours tor inviscai radial diffuser model. 





figure 7. - Multigrid convergence histories for turbulent centrifugal 
impeller model. 



Figure 8. - Computational grid for 6:1 pressure ratio centrifugal impeller. 
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Figure 9. - Radius (RMS Pi and stream surface thickness (BESP) for 
6:1 centrifugal impeller. 



Figure 10. - Static pressure distribution for 6.-1 centrifugal impeller. 








Figure 13. - Computational grid for 6:1 pressure ratio radial diffuser. 
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Fioure 14. - Radius (RMSP) and stream surface thickness (BESP) for 6*1 
radial diffuser. . 




Figure 15. - Static pressure distribution for 6:1 radial diffuser. 
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Figure 17. - Mach number contours for turbulent flow in 6:1 radial 
diffuser. 
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